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Abstract 

We consider the interplay of vaccination and migration rates on disease persistence in epidemiological sys- 
tems. We show that short-term and long-term migration can inhibit disease persistence. As a result, we 
show how migration changes how vaccination rates should be chosen to maintain herd immunity. In a sys- 
tem of coupled SIR models, we analyze how disease eradication depends explicitly on vaccine distribution 
and migration connectivity. The analysis suggests potentially novel vaccination policies that underscore the 
importance of optimal placement of finite resources. 
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1. Introduction 

Countries are increasingly connected by travel and economics. Due to economic disparities and political 
turmoil, extreme heterogeneity exists in childhood vaccination coverage across the two sides of multiple 
national boundaries. It has been suggested that the immunization coverage of neighboring countries or those 
countries well connected by travel can or should be used when crafting national level immunization policy. 
In the case of hepatitis B, Gay and Edmunds yj] argue that it would be four times more cost effective for 
the United Kingdom to sponsor a vaccination program in Bangladesh than to introduce its own universal 
program. When indigenous wild polio virus was eradicated in all but four endemic countries in 2005: India, 
Nigeria, Pakistan and Afghanistan, it was exported from northern Nigeria and northern India and subse- 
quently caused > 50 outbreaks and paralyzed > 1500 children in previously polio-free countries across Asia 
and Africa Q. And in 2007, the WHO estimated that there were 197,000 measles deaths, despite the 82% 
worldwide vaccination coverage. In countries where measles has been largely eliminated, cases imported 
from other countries remain an important source of infection [3]. It is clear that a country needs to be 
concerned with the vaccination rate of a neighboring country as well as its own. 

On another scale, vaccination policies must also take into consideration the subpopulation dynamics 
within a country. Wilson, et al. J3l models linked urban and rural epidemics of HIV and discusses how to 
optimize a limited treatment supply to minimize new infections. Cummings et al. uses data to identify 
a distinct pattern in the periodicity of measles outbreaks in Cameroon before the widespread vaccination 
efforts of the Measles Initiative. The southern part of Cameroon experienced a significant measles epidemic 
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approximately every three years. In contrast, the three northern provinces contend with annual measles 
epidemics. In 2000 and 200L these cyclic outbreaks coincided, exacerbating the situation and causing a 
much more severe epidemic Jj§]. 

Noting that a small contribution of infections from one population to another could drive a new type 
of epidemic that would not normally occur, we study how migration between populations could change 
dynamics and respective herd immunity levels in metapopulation models. We analyze a model of a disease 
imported between subpopulations of a region by short-term and long-term migration with limited vaccination 
coverage. Our initial study is based on the analysis of a system of canonical SIR compartmental models. 
The system allows the rigorous proof of the qualitative affects migration has on herd immunity. The model 
can be enhanced to include more compartments or seasonal forcing, but most of these systems will require 
numerical exploration of trends in spatial synchrony and bifurcation analysis, which will be explored in 
future papers. In this article, we revisit the fundamental ways migration is modeled in metapopulation 
models and how it fundamentally affects herd immunity. 

Migration is often treated as a phenomenological input to maintain incidence in a population that might 
experience local fade-out Jf^. Long-term migration has been analyzed by Liebovitch and Schwartz Q], 
with a thorough derivation of the linear flux term coupling the patches. This approach also agrees with the 
classes of models proposed by Sattenspiel and Dietz and Lloyd and Jansen @]. Keeling and Rohani fioll 
investigated the spatial coupling of dynamics exhibited in models using multiple formulations of migration 
including mass-action coupling and linear flux terms. However, they did not explore the impact of coupling 
in the presence of vaccination needed to maintain disease free states. Additional analysis of mixed long-term 
and short-term migration in transport-related disease spread can be found in 111 llll2lll3ll . These papers derive 
the global asymptotic stability of the disease free state for a new disease. Because there is no vaccination, 
the papers conclude that it is essential to strengthen restrictions of passenger travel as soon as the infectious 
diseases appear. 

Our paper considers how migration directly affects the vaccination levels needed for herd immunity 
against a known disease and how that would impact optimum usage of limited vaccination supplies. We 
investigate the dynamics of models that include mass-action coupling, an assumption that assumes mixing 
occurs at fast time scales, and linear migration, which is more consistent with mixing occurring at long 
time scales. The organization of this paper is as follows: We introduce a coupled compartmental model in 
Section|2]and perform stability analysis of the disease free state as a function of the migration and vaccination 
rates. We also consider normal forms of the bifurcations created by the short-term and long-term migration 
dynamics. Section [3] describes how vaccination rates should be adjusted with respect to short-term and 
long-term migration levels to preserve herd immunity. Section |4] has a summary of our observations and 
conclusions. 



2. The Model 

We start with the classic Susceptible, Infected, Recovered (SIR) model. Let S, I, and R denote the 
number of people in each of the disease classes for a population of size N. Let the parameters /? > denote 
the contact rate, /i > denote the birth/death rate, and k > denote the recovery rate. The vaccination rate, 
< v < 1, represents the removal of a percentage of the incoming newborn population to recovered. The 
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standard form for this system is 



dS PSI 

"4- = ^ -«/-„/, (1) 
dt N F w 



dt 



V/J.N + kI - fiR. 



The death rates in the classes balance the births so that the population size N > is constant. For a detailed 
analysis of the single patch formulation of this system, see Hethcote 1 1-411 . 

We now consider two coupled subpopulations where the disease dynamics of each population are de- 
scribed by the SIR model. Let S k, h, and Rk denote the number of people in each of the disease classes, 
> denote the birth/death rates, and Vk denote the vaccination rates of subpopulations Nt for k — 1, 2. To 
model long-term movement (linear mixing), let c\ > denote the rate of migration from population two to 
population one and vice versa for the rate c 2 > 0. To model short-term movement (mass action mixing), let 
< C3 < 1 be a scaling of the number of infectives from one population who move into the other population 
for a short time and mix with the susceptibles to produce additional infections. Because /3 is proportional to 
the average number of contacts a person can make per unit time, we distribute the contacts for the suscep- 
tibles between the infected people by mass action within and outside the population by using the prefactors 
C3 and (1 - C3) respectively as in Keeling and Rohani lioh . The coupled two population model is as follows: 

— = (l-VilfrNi /iiSi +C1S2-C2S1, 

dt Nt N\ 

dh (l-c 3 )^5i/i ci/3S i/ 2 . 

— = + — kI\ - Hih +c\h -c 2 h, 

dt N\ Ni 

dRi 



dl 



v 1 iu 1 N 1 +kI 1 -iuiR 1 +c 1 R 2 -c 2 Ri, (2) 



dS 2 (l-c,)/3S 2 I 2 

— - = (l-v 2 )fi 2 N 2 /i 2 S 2 + c 2 Si - ci5 2 , 

dt N 2 N 2 

dl 2 (l-c 3 )/3S 2 I 2 c 3y 85 2 /i . 

— = + — kI 2 -pi 2 l 2 + c 2 I x -c\h, 

dt N 2 N 2 

dRo 

= v 2 )x 2 N 2 + kI 2 — )x 2 R 2 + c 2 R\ — c\R 2 . 

dt 

We keep the number of people in the subpopulations constant by letting p = N 2 /N\ and setting the 
constraint c 2 = c\p. This system is overdetermined by the subpopulation constraints, S k + h + Rk — Nt for 
k — 1,2, and therefore the analysis omits the variables Rk for k — 1,2. 

Motivated by the distinct subpopulation dynamics of Cameroon described in Cummings et al. fl, nu- 
merical simulations will use parameters based on Cameroon demographics. The values are listed in TableQ] 
The subpopulation sizes are totals for the northern and southern regions based on data in 15|]. The birth/death 
rates are averages over the northern and southern regions based on data in J^. The recovery rate is a pa- 
rameter that is derived from the biological characteristics of measles. The contact rate was estimated using 
the average age of incident measles cases over the period 1998-2006 from passive surveillance data J5j]. The 
specific results here are fairly insensitive to changes to /?. An SIR model is used here without the exposed 
class but we expect the inclusion of an exposed class would not substantively change our qualitative results. 
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Table 1: Parameter Values for Model Based on Cameroon Data 



Parameter 


Value 


Unit 


Description 


Ni 


4,451,000 


people 


Northern subpopulation size 


N 2 


10,212,000 


people 


Southern subpopulation size 


P 


2.2943 


none 


Ratio ofN 2 /Ni 


P 


700 


year -1 


Contact rate 


K 


100 


year -1 


Measles recovery rate 


Ml 


.0428 


year" 1 


Birth and death rate for Ni 




.0329 


year -1 


Birth and death rate for N 2 



2.1. General system analysis 

We start with a general analysis of the system to determine the conditions necessary for the populations 
to be disease free. 

Proposition 1. System (0 has a disease free equilibrium (DFE) and is given by 

(.S 1 ,I u S 2 ,l2) = (NiSuO,N 2 $ 2 ,0), (3) 

for 

A (1 - Vl)(jUlCi +fXlfl 2 ) + (1 - V 2 )fl2ClP ... 
a l - ; ; , (4) 

n\c\ + +mc\p 

e (1 - Vl)jUlCl + (1 - V 2 )(jlifX 2 + fi 2 Cip) 

J2= ■ (5) 

H\c\ + fiifi 2 + p 2 c\p 

Note that the DFE does not depend on the short-term migration parameter C3. Without long-term migra- 
tion (ci = 0), the DFE simplifies to (Si,Ii,S 2 ,I 2 ) = (N\(l - Vi), 0,^(1 - v 2 ),0), which is the steady state 
for the uncoupled system. Also, there is no steady state for which the disease dies out in only one of the two 
subpopulations if ci > 0. 

The local stability of the DFE can be determined by the eigenvalues of the Jacobian of the system 
evaluated at the DFE. The resulting eigenvalues are 

A x = +cip+n\ +fi 2 + yjffil +C\p+Hx - n 2 ) 2 -4ci(^i -yU 2 )j, (6) 

X 2 = +c\p+fii +fx 2 - yj(ci +cip+pi -/j 2 ) 2 -4ci(yUi -p 2 )j, (7) 

^3 = -i(c 1 +c 1 p+ A / I + / U2+2^-(l-C3)yS(5 I +5 2 ) + VW), 



(8) 



= -i(c 1 +c 1 p+^ I +^ 2 +2^-(l-c 3 )yS(5 I +5 2 ) - Vw), (9) 



for 



W = 4f3c 3Cl (s 2 + p§ !) + 4(p 2 cjS iS 2 + c 2 p) + ((I -c 3 )j3(§i - S 2 ) + a - c x p - m +(i 2 ) . (10) 
The DFE is locally stable if the maximum value of the real parts of this set of eigenvalues is negative. 
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Figure 1: Contour plot of A4 values as a function of ci and C3. Parameters are given by the values in Table [T] with vi = 0.82 and 
i>2 = 0.90. The DFE is unstable for A4 > 0, which occurs for smaller values of c\ and C3. 

Proposition 2. The eigenvalue A4 determines the local stability of the DFE. 
Proof. We can show A\ and A.% are always negative. First, let 

0i = c\ +cip + fii +fi 2 , 

(ID 

62 = (ci +cip + iii -fix) z -4ci(//i -m). 

If we consider 62 as a quadratic expression in c\ with leading coefficient (p + l) 2 , it attains an absolute 
minimum at d02/dci = 0. Solving this equation gives c\ — (pi - 112) lip + l) 2 - Substituting this expression 
into 62 to find the absolute minimum gives 4p(yUi - p2) 2 /(p + l) 2 > 0. Therefore 62 > for all c\, which 
implies A\ and A2 are real valued. 

Upon inspection we see that 6\ > 0, 0\ + y/fh > 0, and therefore Ai < 0. For A2 < 0, we require 
61 > yf&2- This is equivalent to c\ > -/Ji^/ipi + M2P)- Since we assume c\ > 0, this is always true. 
Therefore 6\ - V#2 > 0, which implies A2 < 0. 

For our parameter assumptions, we see that W > by inspection. This implies A3 and A\ are real valued 
and A4 > A3. The only way for the DFE to be unstable is for A3 > or A4 > 0. Because of the ordering, a 
sign change would have to happen for A4 first. Therefore, A4 determines the stability of the DFE. □ 

To quantify how each migration type effects the stability of the DFE, we can monitor the sign of A4 as we 
vary c\ and C3. As an example, we show a contour plot of A4 in Fig.Q]using the parameters in Table Q] with 
vi = 0.82 and V2 = 0.90. When A4 > 0, the DFE is unstable. From the figure, you can see that as c\ and C3 
decrease, A4 increases and the DFE becomes unstable. We now explore the underlying conditions necessary 
in each subpopulation for which migration can cause the die out or invasion of a disease. 

In the absence of migration, we recover the basic reproductive numbers scaled by vaccination for each 
subpopulation for the uncoupled system lfl4}l . Specifically when c\ = C3 = 0, 

6 - vi) fi fli - va) .... 

«i(vi) = and R 2 (v 2 ) = (12) 

K + jl\ K + fl2 



5 



We omit the arguments for for k — 1 , 2 from here on, unless otherwise specified. The basic reproduc- 
tive number is the quantity that defines the threshold between disease absence and persistence, and for the 
canonical SIR model without vaccination Ro = We write these expressions as a function of the vacci- 
nation rate in the subpopulation noting that the inequalities /?# < 1 for k — 1 , 2 implies that the DFE in each 
subpopulation is locally stable. At R^ = 1 for k = 1,2, these two expressions also represent transcritical 
bifurcations for the uncoupled system. As an example, for the parameters used in Fig. Q] /?i(0.82) > 1 and 
the disease would be endemic in N\. Conversely, /?2(0.90) < 1 and the disease would die out in N 2 . This 
motivates us to examine the effect migration has on a simple system with a transcritical bifurcation in each 
component. 



2.2. Normal Form for Linear Mixing 

To understand how long-term migration directly affects the stability of the DFE, we rewrite the system 
as the normal form of a transcritical bifurcation with linear coupling. Since the SIS model has the same 



topology near the DFE as the SIR model, we consider the standard SIS model with births and deaths II15H 



ds 

— = n — Bsi + Ki — us, (13a) 
dt 

di 

— = [isi - id — ui, (13b) 
dt 

with nondimensional variables representing percentages of the population. Since s + i = 1, the system is 
overdetermined and we need only to solve di/dt. 

In the first subpopulation, let x = i and 1 - x = s. Substitute these variables into Eq. (II 3bb and rescale 
time by j3. We repeat the process using the variable y to represent the second population. By adding the 
linear migration terms to these equations, we find 

x — r\x — x 2 — ax + ay, (14) 
y — fiy — y~ — ay + ax. 

Here, the bifurcation parameter r^ = (jRj - 1)/a* for k = 1,2 from Eq. (1121 1. In addition, we rescaled the 
long-term migration rate as the parameter a = c\ //3. 

The steady state (x,y) = (0,0) is equivalent to the DFE in the full system in Eq. ©. In the absence 
of coupling (a = 0), a transcritical bifurcation will occur in the x system, transferring the stability from 
x = to x = r\ at r\ — 0. The dynamics are similar for y, respectively. Linearizing about the steady state 
(x, y) = (0, 0) yields two eigenvalues, 



A, = - (n + r 2 - 2a + ^Aa 2 + (n - r 2 ) 2 ) , (15) 
= - + r 2 - 2a - ^4-a 2 + {r\ - r 2 ) 2 ) . 



2 

A2 



The following analysis uses this linearization approach to conclude when long-term migration can change 
the stability of this steady state. 

We start by considering the case of two isolated endemic populations, which have basic reproduction 
numbers greater than one. We ask if it is possible to stabilize the die out state through the coupling parameter, 
a. 

Proposition3. lfr\,r 2 > 0, then the fixed point (x,y) — (0,0) is unstable for all a e [0,oo). 
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Proof. Upon inspection, Ai is the dominant eigenvalue. Assuming r\,r 2 > 0, then Ai > implies 

Ol + r 2 ) + ( V4ff 2 + (n - r 2 ) 2 ) > 2a. (16) 
Squaring both sides and simplifying, we find 

r 2 + r\ + Oi + r 2 ) ^4a 2 + {r x - r 2 ) 2 > 0, (17) 
which is always true. Therefore, (x,y) = (0, 0) is unstable for all a e [0, oo). □ 

We can interpret this abstract result in the original system by concluding that for two isolated endemic 
populations, the amount of long-term migration is irrelevant to the persistence of the disease. The stability 
of the DFE cannot be changed by migration and intervention by vaccination is necessary for disease die out. 

Next, consider the case where we have sufficient vaccination so that one of the basic reproductive num- 
bers is less than one, while the other is not. Again, we ask under what conditions the coupling parameter 
can stabilize the die out state. 

Proposition 4. Without loss of generality, we assume r\ > and r 2 < 0. Case 1: If —r 2 < r\, then the fixed 
point (x,y) — (0, 0) is unstable for all a e [0, oo). Case 2: If —r 2 > r\, then there exists some a* e [0, oo) 
such that the fixed point (0, 0) is stable for all a > a*. 

Proof. In both cases, assume r\ > and r 2 < 0. 

Case 1: For —r 2 < r\, A\ > reduces to the relationship in Eq. ([T7V This is always true for < r\ + r 2 
and we conclude (x,y) = (0, 0) is unstable for all a e [0, oo). 

Case 2: For -r 2 > n, A] < reduces to the relationship 

( V4a 2 + (n - r 2 ) 2 ) <2a- (n + r 2 ). (18) 
Because both sides are positive, we can square both sides to find 

a > > 0. (19) 

r\ + r 2 

There exists an a* = r\r 2 j(r\ + r 2 ) for a* 6 [0, oo). Therefore, for a > a*, it follows that Ai < and 

(x,y) = (0,0) is stable. □ 

This result implies that in a system with one population supporting an endemic state, there is a minimum 
amount of migration necessary for the system to achieve stability of the DFE. In fact, we can interpret the 
requirement of -r 2 > r\ as y — in the uncoupled system is more stable than x — 0. Therefore, y is sharing 
its extra stability with x. 

For completeness, we can also show that long-term migration cannot change the stability of a stable die 
out state for two isolated populations that have basic reproduction numbers less than one. 

Proposition 5. lfr\, r 2 < 0, then the fixed point (x,y) — (0, 0) is stable for all a € [0, oo). 

Proof. For n, r 2 < 0, Ai < reduces to Eq. ( fT8l . Because both sides are positive, we can square both sides 
to find 

a>^. (20) 
n + r 2 

Since j^- < 0, Ai < for all a e [0, oo) and (x,y) = (0, 0) is stable. □ 

We conclude that long-term migration has a positive effect on the stability of the DFE. The mixing in 
all classes diffuses the force of infection, making it harder for the disease to persist. In applications where 
migration is common, this effect might be significant. 
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2.3. Normal Form for Mass Action Mixing 

To capture the effect of short-term migration for each subpopulation in Eq. ©, we follow the construction 
of model for linear mixing by substituting x and y for i in Eq. (|13bt . In this system, we use mass action 
coupling with the parameter cr - controlling the mass action mixing strength. Specifically, the term 
ail — x)y represents the infectious person from y coming into contact with a susceptible from x. The system 
take the form 

. = (^L^.^d.^ (21) 

V 1 - cr ) 1 — cr 

(r 2 -cr\ 2 0" 

y = h ]y-y — (i-y)x. 

\ 1 - cr / 1 - cr 

Again, the bifurcation parameter rj - (R^ - l)/Rk for k — 1,2 from Eq. (fT2l and time has been rescaled by 
y6(l-cr). 

Performing the linearization about the steady state (x, y) = (0, 0) yields two eigenvalues, 

A i = ^ : (n+r 2 -2cr+ V4cr2 + (n-r 2 ) 2 ), (22) 
2( 1 - cr) v ' 

A 2 = 1 A n +r 2 -2cr- V4cr^ + (r 1 - r 2 )A . (23) 

2( 1 - cr) v ' 

The eigenvalues for this system are of a similar form as those for linear migration in Eq. ([TBI , but multiplied 
by 1/(1 - cr). Therefore, the mass action mixing can change the stability of the DFE as a function of basic 
reproduction numbers in the same settings as linear mixing. 



3. Vaccination responses 

This section directly considers how the migration rates change the vaccination levels necessary to keep 
the DFE stable, which implies the occurrence of herd immunity. It explores whether neglecting the short- 
and long-term migration rates overestimates or underestimates the minimum vaccination rates necessary for 
disease fade-out. 

We first restrict our attention to long-term migration only; i.e., letting C3 = in Eq. (0. For c\ > 0, this 
is equivalent to identifying the bifurcation points in (vi , v 2 ) when A4 — in Eq. ((9). Notice that if c\ =0, the 
subpopulations are isolated. The vaccination levels needed in each for the disease to die out is equivalent to 
solving for Vj and v 2 in R\ 2 < 1 from Eq. ( TT2l . The constant solutions for R^2 = 1 using parameter values 
in TableQ]are shown in Figure|2]as solid black lines, and the disease will die out in both populations in the 
top right quadrant. 

As we increase c\, the boundary for the region of stability for the DFE spreads away from the c\ — 
case, increasing the die out region. The limit is a line, shown by the dotted black line in Figure [2] 

Proposition 6. As c\ — * 00, the bifurcation curve bounding the stable region approaches the line 

, (Mi +H2P)(P-k) (Mi +V2P) 2 
v 2 = — I vi H . (24) 

\H2pj PH2P pH2p(l+p) 

This line is decreasing in v\ , with the slope depending on a ratio of the birth rates and subpopulation sizes. 
Note that if we use the basic reproductive numbers for the isolated subpopulations, Eq. (l24l is equivalent to 



-I — |v 

\P2P 



(25) 

(1+P) \\ fl,(0) ) \ R 2 (0) ) J 
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Figure 2: Boundary of the region of stability for the DFE as we vary ci . The curves represent A4 = using the parameter values 
in Table [T] and C3 = O.The vertical and horizontal black lines represent the vaccination rates necessary for a stable DFE in isolated 
populations (c'i = 0). The dashed line represents the limiting curve as we increase C[ . 



As ^i(O) and/or ^2(0) increase, the V2 intercept increases and shifts the line up vertically. Therefore, the 
attainable stable DFE region in (vi, V2) space decreases, as expected. 

Next, we consider only short-term migration, i.e. letting c\ = 0. Similarly, for < C3 < 1, this is 
equivalent to identifying the bifurcation points in (vi, V2) when A4 = 0. We fix the parameters to the values 
in Table Q] and vary C3 to see the changes to the boundary of the stable region. The solution for C3 = is 
shown as solid black lines in Figure [3] Again the disease will die out in both populations in the top right 
quadrant. As we increase C3, the boundary for the region of stability for the DFE smoothly pulls away from 
top right quadrant, increasing the die out region. 

Proposition 7. The limit as we increase C3 — > 1 is 

_ , _ (K + /U2) +^1) 

V2 ~ (1 - vi)P 2 ' 

Therefore, in both cases, underestimating the migration between populations causes an overestimation 
of the vaccination levels needed for herd immunity. 



4. Conclusions 

In this paper, we consider the effects of short- and long-term migration in coupled population models in 
the presence of vaccination. We study the interplay between the independent vaccination and migration rates 
across different populations. We conclude that neglecting migration effects overestimates the vaccination 
levels necessary to achieve herd immunity. 

We have proven that if two isolated populations support an endemic state simultaneously, migration 
cannot change the stability of those endemic states. Analogously, this is also true for two populations with 
stable disease free equilibria. In contrast, migration can lead to disease die out in the mixed case. If a single 
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Figure 3: Boundary of the region of stability for the DFE as we vary C3. The curves represent A4 = using the parameter values 
in Table [T]and c\ = 0. The vertical and horizontal black lines represent the vaccination rates necessary for a stable DFE in isolated 
populations (C3 = 0). 



population has a vaccination rate sufficient for herd immunity in isolation, low levels of migration from a 
population that is endemic will not necessarily make the disease endemic in both. In fact, increased levels 
of migration can lead to disease die out in both populations. However, migration rates are only physically 
realistic when they are small. 

Our results suggest more efficient vaccination strategies may be identified for groups of countries with 
significant migration between them. For example, instead of increasing the vaccination levels in a population 
that has already achieved herd immunity, sending vaccine to the less vaccinated neighboring country could 
have a greater impact on outbreak levels. The most efficient control algorithm would be to target the stable 
die out region as shown in Figure |2]or Figure [3] 

Conversely, populations for which vaccine delivery is difficult may benefit to a degree by vaccination of 
neighboring countries. More specifically, consider decreasing the migration rates to a country with a lower 
vaccination rate. We show in Figure|4]a policy where N2, which has a vaccination rate V2 = 0.9, decreases 
the long-term migration rate c\ with N\, which has a vaccination rate vi = 0.7. The short term migration 
rate is held constant at C3 =0.1. The decrease in number of new infections for N2 is a small percentage of 
the increase in infections in N\, and we conclude that the policy meant to help N2 has unintended negative 
consequences for N\ . 

In future directions, the model can be extended to include the effects of seasonality. A similar analysis of 
stable periodic behavior can reveal the sensitivity of synchronization to short-term and long-term migration. 
For example, the work of Schwartz lfl6l predicts new periodic orbits that can be excited by the mass action 
coupling in models with seasonal forcing. Specifically, these orbits exhibit long period outbreaks in small 
populations due to mass action coupling. When applying time dependent vaccination schedules, other pa- 
rameters must be considered in addition to the average vaccination rates, such as pulse frequency and phase 
with respect to periodic application. Other techniques can be extended to migration models with vaccine 
control, such as prediction of future outbreaks as reported in Schwartz, et al. 11711 . 
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Figure 4: Time series of infectives in both populations using the parameter values in Table[T] with V\ = 0.7, V2 = 0.9, and C3 = 0.1. 
The value of c\ is decreased to the constant noted in each window. 



Acknowledgments 

We gratefully acknowledge support from the Office of Naval Research. The authors were also supported 
by the National Institute of General Medical Sciences (Award No. R01GM090204). DATC holds a Career 
Award at the Scientific Interface from the Burroughs Wellcome Fund and received funding from the Bill 
and Melinda Gates Foundation Vaccine Modeling Initiative. The content is solely the responsibility of the 
authors and does not necessarily represent the official views of the National Institute of General Medical 
Sciences or the National Institutes of Health. We would also like to thank Leah Shaw and Luis Mier-Y- 
Teran for their useful discussions. 



References 

[1] N. J. Gay, W. J. Edmunds, Developed countries could pay for hepatitis B vaccination in developing 
countries, British Medical Journal 316 (1998) 7142. 

[2] B. Aylward, R. Tangermann, The global polio eradication initiative: Lessons learned and prospects for 
success, Vaccine Suppl 4 (201 1) D80-5. 

[3] World Health Organization (WHO), Measles Fact Sheet No. 286, December 2009, 
http : //www . who . int/mediacentre/f actsheets/f s286/en// 

[4] D. P. Wilson, J. Kahn, S. M. Blower, Predicting the epidemiological impact of antiretroviral allocation 
strategies in KwaZulu-Natal: The effect of the urban-rural divide, PNAS 103 (2006) 14228-14233. 

[5] D. A. T. Cummings, W. J. Moss, K. Long, C. S. Wiysonge, T. J. Muluh, B. Kollo, E. Nomo, N. D. 
Wolfe, D. S. Burke, Improved measles surveillance in Cameroon reveals two major dynamic patterns 
of incidence, International Journal of Infectious Diseases 10 (2) (2006) 148-155. 



11 



[6] N. M. Ferguson, C. A. Donnelly, R. M. Anderson, Transmission dynamics and epidemiology of 
dengue: Insights from age-stratified sero-prevalence surveys, Phil. Trans. R. Soc. London, Ser. B 354 
(1999) 757-768. 

[7] L. S. Liebovitch, I. B. Schwartz, Migration induced epidemics: Dynamics of flux-based multipatch 
models, Physics Letters A 332 (2004) 256-267. 

[8] L. Sattenspiel, K. Dietz, A structured epidemic model incorporating geographic mobility among re- 
gions, Mathematical Biosciences 128 (1995) 71-91. 

[9] A. L. Lloyd, V. A. A. Jansen, Spatiotemporal dynamics of epidemics: Synchrony in metapopulation 
models, Mathematical Biosciences 188 (2004) 1-16. 

[10] M. J. Keeling, P. Rohani, Estimating spatial coupling in epidemiological systems: A mechanistic ap- 
proach, Ecology Letters 5 (1) (2002) 20-29. 

[11] J. Cui, Y. Takeuchi, Y. Saito, Spreading disease with transport-related infection, Journal of Theoretical 
Biology 239 (3) (2006) 376 - 390. 

[12] Y. Takeuchi, X. Liu, J. Cui, Global dynamics of SIS models with transport-related infection, Journal of 
Mathematical Analysis and Applications 329 (2) (2007) 1460 - 1471. 

[13] J. Liu, Y. Zhou, Global stability of an SIRS epidemic model with transport-related infection, Chaos, 
Solitons & Fractals 40 (1) (2009) 145 - 158. 

[14] H. W. Hethcote, The mathematics of infectious diseases, SIAM Review 42 (2000) 599-653. 

[15] J. D. Murray, Mathematical Biology, Springer, Berlin, 1989. 

[16] I. B. Schwartz, Small amplitude, long period outbreaks in seasonally driven epidemics, J. Math. Biol- 
ogy 30(1992) 473-491. 

[17] I. B. Schwartz, L. Billings, E. M. Bollt, Dynamical epidemic suppression using stochastic prediction 
and control, Physical Review E 70 (2004) 046220. 



12 



